m 



Spectrally determined singularities in a potential with an inverse square 

initial term 

Deinetrios A. Pliakis 
Department of Electronics, TEI Crete 
Romanou 3, Chalepa, Chania, GR-73133 
O' pliakis@chania.teicrete.gr 

(N 
;-H ■ April 10, 2013 

< 

Q> . Abstract 

We study the inverse spectral problem for Bessel type operators with potential v{x): H,^ = —d^ + 
(~| ' -^ + v{x). The potential is assumed smooth in {0,R) and with an asymptotic expansion in powers 

QhI and logarithms as x ^ 0^,v{x) = 0{x°'),a > —2. Specifically we show that the coefficients of the 

fH ' asymptotic expansion of the potential are spectrally determined. This is achieved by computing the 

expansion of the trace of the resolvent of this operator which is spectrally determined and elaborating 
the relation of the expansion of the resolvent with that of the potential, through the singular asymptotics 
lemma. 



C^ 



> . 1 Introduction 

oo 
oo 
^o , 

fS| , Let R > 0,a > —2 and t; be a real valued function, smooth in C^((0,i?)) that has an asymptotic 
.^ ' expansion together with all its derivatives as x — ?• 0+ of the form 

o' 

^^ ■ v{x) ~ x"^"' -''" 



X 



,VnX 
n=0 



where {kn} is an increasing sequence of positive real numbers. This function appears in the singular 
Sturm-Liouville operator H = —d^ + -% + v{x) on (0, R) with Dirichlet boundary conditions at 0, R. In 
^.' this paper we examine whether the spectrum of this Sturm-Liouville operator uniquely determines the 
sequence of asymptotic coefficients {vn} of the potential function. 

We prove that this is indeed the case provided that i^o 7^ 0; this is achieved by the explicit construction 
of the spectral invariants that in turn determine the asymptotic coefficients of the potential function v 
from the spectrum. The function 

is defined for small 77 > and is of course spectrally determined. The 77—7-0 asymptotic expansion of / will 
provide the desired spectral invariants: the singularity of the Sturm -Liouville operator H manifests itself 
in the appearance of coefficients in the 77 — )• 0^ asymptotic expansion of fHiv) that depend polynomially 
on the coefficients {vn} of v. Precisely we have that this asymptotic expansion is of the form 

00 00 00 

/H(r/)~ J]A„772"+i + 774( J] i?„772"+Mog77 + j;C„77'") (1) 

n=0 n=— 1 n=0 



where In is real but not an odd integer. Evidently it contains apart from pure odd powers also new terms, 
odd powers with logarithms as well as pure powers other than the odd ones. The sequences of asymptotic 
coefficients i?„, C„ of the asymptotic expansion of / are polynomials in the asymptotic coefficients of the 
potential {vn}', it is exactly this fact that allows us to determine recursively these {vn} from the "heat 
invariants" that are equivalent to the r/ — )• 0"*" asymptotic coefficients of fHiv)- The exact formulation of 
this tehorem is the following: 

Theorem. Let x be a smooth cutoff function, x = 1 ^'^ ^ neighborhood of the origin, then the function 
fxtniv) ^^-5 (^^ asymptotic expansion of the form 

oo oo oo oo 

fxMv) ~ E ^^k+iv'"^' + E ^-^'" + E C2fc^'' + E ^2fe+ir/2'+' log V 

fc=0 n=0 fc=0 fc=0 

where z^ £ R \ '^,Bk,Ck,Dk are homogeneous polynomials in the asympotic coefficients {vn} of the 
potential of degree z if we define deg{vn) = kn + 2. 

Inverse spectral results of the form we described above appeared first in [4], [5], [6], [10], [9]. Actually in 
[9] certain coefficients are calculated in the course of calculations of operator determinats and in the rest 
the generic case is studied exhaustively. The general case requires the treatment of certain exceptional 
cases. This treatment is provided here and leads to the full result. 

The conical singularities have been treated extensively in the literature for diverse purposes giving 
rise to diverse results and calculations. We refer to the first papers that dealt with asymptotics of either 
the heat [12] or the wave kernel [13] on spaces with such singularities. 

Possible applications of these results are indicated already in [6]. We recall these briefly here. Po- 
tentials of this type arise in the wave equation for a vibrating rod of variable cross section, when the 
cross-sectional area of the rod vanishes quartically (as a function of the distance from the end of the 
rod) at one point. In the same spirit we could determine the profile of a surface of revolution with a 
conical singularity on the axis, asymptotically to all orders, from the spectrum of Laplacian restricted to 
functions with polar symmetry. 

We indicate here briefly two applications in physics. Furthermore the present study exhausts the 
determination of a confining potential which is hydrogen atom-like at short distance [15] , from the complete 
set of bound state energies. This is indicated in [5] but only for the case of the s— wave. Hence the present 
study fills in the higher angular momemtum instances for potentials of asymptotic expansion of the form 
given above. We will present here briefly the construction in [1]. Confinement is expressed by Dirichlet 
boundary conditions on the surface of a sphere. The Hamiltonian with a radial potential is 

H = -A-i + v{\x\) 

where v € r°°((0,i?)) with v{x) = 0{x°'),a > —2 as x — t- 0. The boundary condition on the domain of 
H is tp{R) = 0. The /— th spherical harmonic Hamiltonian is unitarily equivalent to 

Ho = -d'^ + ^-^^ + vir) 

which is an an unbounded operator on L^((0, i?); dr). The boundary conditions inherited from the 3- 
dimensional problem are /(O) = f{R) = 0. A physically meaningful question is whether the potential v 
could be determined from the mass spectrum of the bound states associated to a given spherical harmonic 
because both the mass (the energy) and the angular momentum (the order of the spherical harmonic) are 
measurable in the laboratory. 



The study we perform answers this affirmatively when the potential is a real analytic function on 
(0, R) and that near is written as a convergent series of the following form 

v{x) = x~^ 2_, Vkjx"'' log-^ X, 

for {afc} an increasing sequence of positive real numbers. If the function v is not of this form then we can 
only obtain the full asymptotic expansion of the potential at small distances. Moreover we can deduce 
several corollaries of the preceding theorem that are interesting in physics. For instance availing ourselves 
of the results of [8] on irregular singularities, we state the following: 

Corollary. Let H = —d^. + v(x) + P{x) he an oprator on R-|_ with Dirichlet boundary conditions. 
Let V he a real valued analytic function, rapidly decreasing at oo, while near it is given by a convergent 
series of the form 

v{x) = x~ y, Vkjx'^'' log-' X 

for {ofc} o-n increasing sequance of positive real numbers. The function P is given and is of the form 

N 
P[x) = ^ Ofcx'^" 
n=0 

for < Po < ■ ■ ■ < Pn and a^ > 0. Then the spectrum of H determines v. 

The preceding result incorporates for instance the case of quantum particle systems trapped in mag- 
netic fields. An immediate application of the preceding corollary refers to the 2-dimensional magnetic 
trap [1], due to a vanishing magnetic field at oo. Precisely, let < e < 1 and (r, 9) be the coordinates in 
R^ and consider the quanutm motion in a magnetic field of the type B = {2 — e)r~'^. The Hamilotnian 
for a particle subject to such a magnetic field as well as to an unknown radial potential v{r) has the form 



H = -dl + \dl + r2(i-^) + 2ir~'de + v{r) 



„2 

and restricting to angular momentum m we obtain an 1-dimensional operator of the above type provided 
that v{r) is of the required form. We notice also the example from [17] of a quantum particle moving 
in a magnetic field with fixed angular momentum in the direction of the field. Similar conclusions could 
be obtained in that case as well. The results in [2] allow us to incorporate the case of a pair of opposite 
charges with fixed relative angular momentum moving in a homogeneous magnetic field and interacting 
through an unknown potential that we wish to determine spectrally when the latter satisfies the preceding 
assumptions. 

Actually we insert a cutoff function x-, suppx C [0, R) in order to deal the operator in the half line. 
The 77—7-0 asymptotic expansion of 

ffj{r,)=r,hr{l + r,'Hr' 

- H is the operator in (0, 00)- coincides mod{r]°^) terms with that of fn. 

The paper is organised as follows: we start with the neccessary facts concerning the operator domain, 
the operator estimates and we construct the resolvent of the unperturbed operator. Next we give the 
existence and the precise form of the asymptotic expansion for the operator function we introduced above. 
We conclude with the calculation of the asymptotic coefficients that provide the recursive relations that 
determine the potential asymptotic coefficients. The solution of these require certain elementary estimates 
that are included. 

Acknowledgements. 1 would like to thank my advisor Professor C. Callias who suggested the 
problem as well as Professor E. Floratos who pointed out the possible appliction to the case of magnetic 
fields. 



2 Operator Domain and Operator estimates 

2.1 Function spaces for singular heat expansions 

We will consider operators of the form H = if^ + v, for H^^ = —9^ + -^ and f is a real valued potential 
function belongs to the space r°°(0, R) of functions that have an asymptotic expansion as a; — t- 0"*" together 
with their derivatives. This space is described as follows that are encountered in the resolvent expansions 
that we are going to deal with. 

The space r°°(R+) consists of the C°°(R+)— functions that have asymptotic expansions as x — )■ 0"^, 
together with all their derivatives, in complex powers of x and integral powers of logx. Precisely, / S 
r°°(R+) if / € C°°(R+) and there exists a 5 : C — > Z+ such that Yli^z<a ^i'^) ^ °^ ^^^ each a G R, and 

97 fix) = E E /-J^"(^' log' ^) + 0(x'^+^"™) 

»2<aZ+9j<5(2) 

for some (5 = (5^ > 0, some fzj G C and all a G R. S is called the asymptotic character of the expansion. 
Define the following differential operators on C°°(0,oo), for given S : C — )■ Z_(_ as described above 

Pz[s]= n {xd^-zf^''\ p:[s]=Pz[s]{xd^-zf^^"^= n {xd^-z'f^^'^ 

^z'<Mz ^z'<Mz 

The following proposition is proved in [7]; though it is elementary it provides criterion that shades light 
on these spaces 

Proposition 1. / G r°°(0,oo) iff f ^ C°°(0,oo) and there is an asymptotic character such that 
Pz[S]f{x) = 0{x'-^) for allzeC and e > 0. 

The space r°°((R_^)^) consists of all the functions / G C°°((R)^) for which there exist 81,82 such that 

{xd,r{ydyrP^^[8i]WA82]f{x,y) = 0(x«^^+^^/^^+^^) 
for all 21,^2 £ C,si,S2 G Z+,(3;,y) in any compact neighborhood of 5(R^) and some 6j depending on 

Zl,Z2. 

Evidently {81,82) is analogously the asymptotic character of / G r°°(R^). If / G r°°(R4.), we let 
Dk,jf{0) denote the coefficient of x'' log-* x in the expansion of / as x — )• 0. 

Now let a > —2, {A;„,} be an increasing sequence of positive real numbers then we have that 8{z) = 1 
for z = kn + a and 8{z) = otherwise, equivalently, the potential function has the following asymptotic 
expansion as x — )• 



v{x) ~ x°'/__. 



00 



VnX 
n=0 

Actually for < x < e and a > — 2 the obvious estimate x" < ex~^ + B^ implies for cp G Cg"(R+) and 
c < 1 that 

\\V^L^<c\\H,<t>h^+bUh^ 

which in view of the Kato-Rellich theorem reduces the domain questions of H to those for Hf^. 
Through the Hardy inequality we conclude that if k > — | then the operator H^ is positive and hence 
possesses at least one self adjoint extension, the Friedrichs' extension. Recall that if k > |, then Weyl's 
criterion implies that the operator iJ^ is essentially self-adjoint and hence the Friedrichs' is the only 
extension. The domain of self adjointness consists of the functions 

V{H) = {^€ l2(r7), 1192^11^2 < 00, IIx^VIIl^ < 00} C l2(R7; ^) n ^^(R^). 



The Sobolev embedding theorem imphes that these are L^-functions with absolutely continuous first 
derivative. Additionally, since (j) G T^{H) then ||-^|| < oo which in turn expresses the Dirichlet boundary 
conditions, (^(0) = Wvux^q (j){x) = 0. 

2.2 The resolvent. 

Let H = —d'^ + f be a Shrodinger operator, with real valued potential function. Clearly, it is formally 
symmetric. If (j), ip are the unique elements of kerH that are integrable at 0, oo respectively then the 
resolvent of H is given by 



w,{4>{X),i^iX)) 



W,{(t>{X),i;{X)) 



where W is the Wronskian of the solutions of {H—X)(p = 0. In the case of the operator Hf^^, 



q2, k I o 



IS 



the solutions are expressed through the confluent hyper geometric functions jl9| . [16] : 

^{x, X) = M/,,^(2x\/^), ^P{x, X) = Wf,,u{2xy/^) 

where the indices are /i = ^= and u = \/ k + \- Their Wronskia 



W 
The confluent hypergeometric functions coincide for a = with the Bessel functions: 



Mo,^(rEV=A) = 2^T{iy + 1)V^L{ 



-A), 



WoA^^/^) 



iKA^^X), 



the resolvent of H^^^ could be obtained from that of i/^ by Neumann series. The latter will be used in 
the operator estimates that follow. Therefore we form the Whittaker Green's function: 



Rx{x,y) = e{x-y) 
+ e(y-x) 



T{u-fi + i)r(z. - /i + l)M^,,{2yV^)W.f,,,i2xV^) 
2^/^r(2i/ + 2)r(z^ + // + ^) 

T{u-fi + ^)T{u -fi + l)Mf,^^{2yV^)W-f,,u{2xV^) 
2V^r(2i/ + 2)r(z^ + /u + ^) 



2.3 Operator estimates. 

We denote by R^ = (A — Hf^)~^ the resolvent of H^. The positivity of H,^ implies that 
1. The operator norm is 

for |A| — 7- oo uniformly in the cone \^X\ > -(KA + e): 

Since R^L^iR^) C V{H^) we choose then G L'^(R^): R^cp = ip- It follows that for k < -e: 



||(A-i7,)V^||i. = |A-A;| 



2||„;,l|2 _L 
L2 + 



H. 



>\{X-k)\' 



2 
L2 



+ 2(3f?(A:-A)(V,^«V)) 



and the estimate follows. 



^ Since our study is effected in the sector of |9A| < ■j-(5RA + e) then we assume that v — ^ > — i 



The resolvent R^ is represented by the kernel for a = -^^'- 

a a a a 

We have the following classical estimates for the classical trace norms ||^||fc := tr{\A\^Y'^: 

2. Let (j) e C^(R) then H-R^^Hi = 0{\\\-^/'^) for A in the cone |9A| > \{^\ + e). 

Proof: For A < — e and v in the above cone it holds that 

\\R'Mi<\\{l + {\-u)Rl)\\ll\\Rl<l^\\i 

as well as that for cr^ = —-r) : 

\\Rl(t)\\i= / dxx(^{x){K,I,){-). 
Jo cr 

Hence we have the required estimate for 

X 

dxx(t){x)[KyI^){—) < a\\x(j)\\L2\\K^I^\\L2. 
o- 



3. Let a G R_ , a + d < 2 then 



\x-d',RlU.=0{\X\'i^-') 



Proof. This follows for d = from the inequality given in the beginning while for d = 2 it results from 
the fact that —d'^< H^^i and therefore 

\\dm < \\h^_,,rW < (ll'^lli. + \x\'\\Rim). 

4 

Furthermore for d = 1 we have for ^ = R^(j) that 

(rE-"c>,.^,x-"a,V') = (x-2-V,-52^) +2a(a,^,x2("-IV') < 
< 2<5||x2-^||i. +2ri52V^||i. +2e||x-(2"+i)^||2, +2e-^||a.V'|li2 

1 a 

the latter choosing 6 = |A|'^~2 and e = |A| a' gives the required estimate. 

2.4 The Neumann series. 

Let a > —2, ko = and {kn}'^=i an increasing sequence of positive real numbers. The potential v has the 
asymptotic expansion at the origin of the form: 



v{x) ~ X°' 2_]vnX^ 



unJ 
n=0 

In the sequel we will distinguish two cases: 

• Case I. Let a = — 1 and vq ^ 0. 

In this case the distributional trace of the resolvent Rx{H), f}j[—i=) := Tr{R\{H)) is computed 
using Neumann series around -ff^.a = — "^x + "^ + f • 

1 ^ 

fH{^=) = y^Tr{Rx{H^^^){v - -)yRx{H^^^) + Tr{Rx{H^,^){v - -))^+'Rx{H). 
V-A n X X 

j=0 



• Case II. Let a > — 1 and vq ^ 0. 

The Neumann series is based on the resolvent of the usual Bessel operator K^. The case of — 2 < 
a < — 1 is treated within the frame of case B. 

The Neumann series around the Bessel operator, due to its behaviour under scaling, is used in order 
to obtain the general form of the asymptotic expansion. 

3 The asymptotic expansion 

In this paragraph we '11 establish the existence and the precise form of the asymptotic expansion of the 
distributional trace of the resolvent Rx{H) of the operator H = —d^ + -^ + v{x) by means of a Neumann 
series based on the resolvent Rx{Hk of H^ = —d"^ + -^, for a cut off function x ^ Co°(R-+))X = 1 in a 
neighbourhood of zero : 

1 ^ 

fxM-7=} ■■= TrURxiH)) = Y,Tr{x{Rx{HM'Rx{H.)) + 0{\\\^^^^) 

i=i 
Denote by Ij[-^=) = Tr{x{R\{HK)vy R\{Hk))- Actually we have the following 

Theorem. Let x ^e a smooth cutoff function, x = 1 in a neighborhood of the origin, then the function 
fx,Hi'n) ^^-5 ^"iT- asymptotic expansion of the form 

oo oo oo oo 

fxMv) ~ E ^2fe+ir?''^+^ + Y. BnV'" + E ^2fc^'' + E D2k+iv''-^' log V 
k=0 n=0 fc=0 fc=0 

where 2^ S R \ Z,Bk,Ck,Dk are homogeneous polynomials in the asympotic coefficients {vn} of the 
potential of degree z if we define deg{vn) = kn + 2. 

The proof will be achieved in two steps: first we prove the existence by appealing to the singular asymp- 
totics lemma that we recall below and in the second step we use a scaling argument in conjuction with 
the asymptotics of the resolvent of the Bessel operator i/^ = —d^ + -^ which as it is well known it is 
represented by the kernel : 

R.{x, y) = ixy)"2[e{x - y)K,{-)h{^) + Q{y - x)h{-)K,{^. 

a a a a 

Step 1. The trace 

Ij =tr{xRx{H^){vRx{H^)y-^) 

is the j + 1-tuple integral 

oo roc roc 



dx j dxi . . . j dxjxix)Rx{x,xi) . . . R\{xj,x) 
Jo Jo 

which is the sum of integrals, each one for an ordering of the the variables x,xi, . . . ,Xj hence for x > 
xi > ■ ■ ■ > Xj: 

/ dx dxi... ~ dx,x\ix)xlv{x^)...x-'^v{x,){K,i-)l4^)fr[K,i^)I,{^) 
Jo Jo Jo (^ (^ fJl cr a 



and we perform the coordinate change (blow up) 

and 

xi = x9i, ... ,Xj = xOj^i, 9k = ti... tk 

that leads to a sum of integrals of the form: 

'^ d9i r^ dd^ ^ - - ^'"^ 

/o ^1 '"io G, ^^^ . „ ^^^ 

which is finally of the form suggested from the SAL 



rdxx'^xix) r^... ['^tlv{xe,)e]K'A-)i'ui-o,)UiK.iu)i-ek) 

Jo Jo »i Jo (^j t-, ' a a H: (^ 



where 



if we set 



and 



r°° dr 
U<^)= -x(x)x3.F(e,x) 

Jo X 

Fitx):= -^■■- ^iB,ic,e,)n,,{C;eu--- ,e,.i)v,ix;eu--- ,ej). 

Jo t7l Jq Uj 

j — l 

j 

•j{x;9i, ■ ■ ■ ,9,) := Il{x9i)\{x9i) = 0(x^'("+2)) 



Now we see that Vj G r°°((R+)) with asymptotic character S + ■ ■ ■ + S since v G r°°(R-|_) has character S, 

j— times 

v{x) = 0{x'^),a > -2. Furthemorethecharacter of (Qi.-P^)(^,6'i,--- ,9j) is S{k) = 1 for A; G J + 1 + 2Z+. 
The estimates for the Bessel functions as ^ — ?• 0: 

From these it follows that F{S,,x) G r°°((R+)2) and also that 

\p:Fic,x)\<ixo'^^+'h,{o 

where we have that 



Rez' <Rez 

Then we conclude that 






I,ia) ^a^f2 ^n'^^'" + ^'""HX; 4'^^'' + X; Cj-'V^'^+I + J2 D['^a'H^og a)) 



n k=0 k=0 k=0 

where the Zn ^ j + 1 + 2Z_|_ and the Bz , Cj^ , D^. are polynomials in the coefficients of the potential of 
the potential of degree k + j if we define degv^ = z + 2. These suffice for the B's in the expansion. 



Step 2. We have to establish the nature of the non classical terms logarithmic terms. This will be 
achieved by induction on the ordered couples (j, k) where j is the order of the term in the Neumann series 
and k is order in the asymptotic expansion of the potential. For k = a the assertion follows from above: 
everything depends on vq\ indeed j = follows from [5]. Assume that k > a. Then consider 

/, = Tr{x{Rx{H^)vyRx{H^) 

Notice that up to 0(A~°°) we could embody the cutoff function in the potentials V. Accordingly let 
A'' E Z, /cat > — 2 and split the potential as v = utv + Rj\[, where vn{x) ~ x°' X]n=o '^nX^"- The following 
commutator identities for A;^ > ^ — 1, r < n and e = kn — kr > give that 

x'^-RxiH^) = x''^Rx{H^)x' + 2ex''^Rx{H^){x'-\xd^ + ^-^)R^{H^), 

xd^Rx{H,) = Rx{H,)xd, - 2Rx{H^fH^ 

and in turn allow us to vary the number of the Rx{H^) factors. This in view of the fact that the asymptotic 
expansion of IuK^{x) as re — )• oo contains only odd powers, in view of \dxRx{Ho) = —2adaRa, leading to 
the polynomial dependence of the coefficients of the logarithms and the even powers: 

oo oo oo oo 

fc=0 n=0 fc=0 fc=0 

where z^i G R \ Z, z > —2 and the B^, Ck, Dk are homogeneous polynomials in the asympotic coefficients 
{vn}^=Q of V of degree z if we define deg{vz) = z + 2. Then summing up we obtain an expansion of the 
same form and we appeal to the results in par. 5 of [1] in order to keep the odd terms in the log's. 
In conclusion we supply the asymptotics of 

as ^ — 7- uniformly in < t < 1. For this we calculate the Mellin transform using the integral represen- 
tation: 

If /-oo 2 J 

cm = KA-^)iA-^) = I ^^^,Mtx)M^)-- 

The Hankel transforms 1201 allows us to obtain further: 



.,. {2tr r x'-+2 f^ j,{xvi+t'^-2ut) dx 

^'^^^^T{u + \)Jo x^ + ^-'Jo (l + t2_2nt)t X- 

From this we get that the Mellin transform is 



V^r(i. + i) ' ' ' Jo ' ' " ' {i + tr 

A Taylor expansion then leads to the series 
where 

h{t)= /\^-^+^'( , ^~2 y'^dii. 

Jo 1 - jT+tyi'n 

9 



This integral satisfies that I^ = 0{k 2 ^) as A; — )• 00, uniformly in t, where (5=^forl>e>0, z^>^ 
and 6 = ly for i^ < 2- Stirling' s formula then suggests that the summand behaves as 0{k~^^ ~ ) and 
hence it is uniformly convergent for < t < 1. Finally we conclude that the asymptotic expansion of 
gt{^) as ^ — >• contains only odd powers of ^. g 

4 The recursions for the coefficients of the expansion 

We will consider the two cases indicated in section 2.4 as follows: 

4.1 Case I. 

The first order term in the Neumann series is reduced mod{X~°^), to 

Setting a = -4= we arrive through the formula 



a2 



at the integral 



where 



hi^r) = —adaIi{(T) 



r°° 1 dr 

h{a) = a. g^,^,(-)Vix)^ (2) 

Jo y ^ 






for M^^y^W^^y being the Whittaker functions y = f ,Ai = ^^^^. This Green's function could be simplified 
in the form that we '11 use in the sequel 



Gi^Av) '■'- 



Qj^Av) 



Gf.Ay) = y'"^'w^Ay)Mt.Ay) 



where 

for the integrals 

a' 

^.,.... 7_,- - -- ^i + .i 

w^,y{y)= / e-^y'\92{i + e2r~H--^)~^d92, 

Jo ^2 

ay{n) = r(i/ - /i + -)r{u + fi + -). 

The study of the a — >■ asymptotics of the above intgeral will provide the inverse spectral result. Indeed 
the Singular Asymptotics lemma provides the asymptotic coefficients of the integral o" — ?■ 0. For brevity 
we introduce the notation 

2 -2 



In order to extend meromorphically it out of the strip —l + iy<^s<uwe differentiate with respect to ^ = jji^ppi- The 
poles as t — >■ 1 remain located there. 

10 



Then we distinguish between the following cases: 

A. 2v G Z-|_, 1/ = n + 2 • For / = 0, . . . we appeal to the logarithmic terms 

21 
j=0 

and also to the even powers that are given 

21 j 

i=o ^' 
Actually for the inverse spectral result we'd like the explicit forms: 

D2i+3,ih{0) = Bli^^{v)v2i+2 -{1 + l)Bli^^{u)voV2i+i - {I + !)(/ + -)Bli_,2{^)viv2i 

+ P2l+2{v) 

and also that 

D2n+2j+2,0h (0) = C^j+2n+l {v)v2n+2j+l + {j + n + -)Clj+2n+l {T^)vQV2n+2j 
+ Q2n+2j+l{v) 

where we have denoted the coefficients 

Bl{n) := DlAv^CjA = Dk-j,oCj,,y, Cl{v) := ul^V^CjA = Uk-jfi{Cj,u) 

and the P, Q are the polynomials suggested precedingly by the general form of the asymptotic expansion. 
Notice that if / = 0, . . . , n then we '11 see in the sequel that: 

therefore coefficients Vk,k = 1, . . . ,2n + 1 are determined immediately. For the remaining we appeal to 
the pair of equations; hence we have to establish, provided that vq ^ 0, that the determinant: 

Afc(l^) = B2n+2k+2{'^)C2n+2k+l{^) ~ B2n+2k+2{^)C2n+2k+l{^) 

does not vanish. 

B. 2i/ G R \ Z_|_. The inverse spectral result requires again the logarithmic odd powers as above, which 
are provided from the first order terms 

D2l+3,llli0) = D%^{U)V21+1 + P2l+2{V) 

as well as the the pure powers a G R, a > —2, a / 0, 2, ... : 

Do^iim = uiA[Po.G^,A-) ■ ^]u-=o- 

In this formula P^ = nz<a(^^^^^^~^) ^^"^ hence we obtain the set of equations 
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Notice that if a 7^ 2zv + 2j + 1 then C^_^2 ¥" ^^^ hence Va is determined by the first terms for a ^ 2z> + A;. 
These in particular contain the coefficients for a = 2z^ + 2j + 1, j = 0, . . . . 



-D2v+2«+4,0^l(0) = C2^+2i+3(^)^0«2!y+2i+2 + C'2i/+2i+3(^)^0^2j/+2«+l + Q2j/+2«+3( 
-D2j.+2«+3,0-^i(0) = C2^_,_2i+2(i^)^2i/+2i+2 + C'2j,+2;+2(l^)^^0«2!/+2i+l + Q2i/+2i+2 



V] 



where Q are the same polynomials that appear also above. Again we have to stablish that determinant: 

Afc(l/) = C2^+2fc+2('^)C2V+2A.-+l('^) - C'L+2fe+2('^)C2i.+2fc+l(^) 

of the coefficients of the preceding system of equations is nonvanishing. 

The Mellin Transforms of Cj^v Finally all these asymptotic coefficients coincide with the coefficients 
of the Laurent expansion of the Mellin transform 

as it is explained in the appendix on the Singular Asymptotics Lemma. In the sequel we' 11 employ the 
identification provided by the identities: 

^fc,o(cj>)(0) = ReSs=-kCj,u Uk,o{'^j,u) = Cj^y{-k) (3) 

The function Cj^yijj) is written under the following change of variables 

^1 = 2^1-1, ^2 = ^if4^ 

1 — O2 

in the form: 

Further we perform the change of variables 

(0, 1) X (0, 1) B [61,62) ^ it ri) G (0, 1) X (0, oo), 

i — [72 

d6id62 _ dCdi] 

1-92 ~ 2r{C,7]){l + 7]) 

where r{^,r]) = (1 — ttx^)^. In these coordinates the preceding integral is written as: 

In the sequel we '11 derive its asymptotics as y — )■ 0. To that end, we perform the change of variable 
^ = YTf to obtain finally by abusing notation for the function r{^,'q) that 

;°° e'^r/^'^+^dr/ f^ ^^'U^ log^' ^ 
CjAy) 



1 + 77 Jo r(e,r/) (1 + eP+i' 
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By a Taylor expansion of the function r(^,r/) we obtain the Melhn transform of the function Cj^^{s): 

oo 

Cj^^{s) = {-iy2'-'^''-^T{-s + 21/ + 1) ^ 4{u)T{s + n + l)T{n - s) 



where 



{2l + l)B{l + ,. + lu + ^] 



d» 



{V{l + l)f 



These combined with Stirling's formula suggests that the series converges absolutely for < Ks < 1 and 
represents a meromorphic function with poles at the integral points. Using the Laurent expansion of the 
Gamma function we obtain that for / < k 

2'-^T{-s + 2u + l)r(s + / + l)r(/ - s) = 
(_l)fc-«r(/ + k + l)T{k + 2z/ + 1) 
2*^-ir(A:-/) 

(^(/ + k + I) + 4){k - I + I) -\og2 + ^(21/ + k + l) + 0{s + k)\ 



s + k 
and for I > k 

2'-^T{-s + 2u + l)T(s + / + l)r(/ -s) = 2-('^+i)r(A; + 2v + l)T{l -k + l)T{l + k + l) + 0{s + k). 

Furthermore let N = [2u] then the Taylor expansions at — (2z/ + k),k € Z+ give that for I < N + k: 

2-^v(-s + 2. + i)r(. + , + 1) . r(/ - ,) = r(;. + 4. + i)r(2. + fc + i + i) ^ 

^ ' ^ / V / sm2z/7rr(2i/ + A:-/) 

and for I < N + k we have 

2'-^r{-s+2iy+i)r{s+i+i)-r{i-s) = 2-'^^''+''+^'>r{k+4iy+i)r{i+i-2u+k)r{2u+k+i)+o{s+2u+k) 

The identities. At this point notice that since 

Mo,.(x) = 2-r(z. + hv^x-'^e^'^Ux), WoAx) = '^—^^^^:^x-''e'K,{x) 

then 

CQ,v{y) = yIu{y)K^{y) 

and this suggests further that 



Co,u{s) 



r(^-| + ^)r(i-§)r(f) 



The general form of the asymptotic expansion suggests the absence of logarithms in the odd powers that 
in turn implies the identities B2^j^^_Au) = 0. In parallel the exact form of Co,;/ suggests that for u = n + i 
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then B^{n + i) = for k = —2{n + I) whereas C'2!y+2fc+3('^) ~ *-*• ^^ '^^ employ later these identities but 

prior to that we study the resulting coefficients. 

The asymptotic coefficients in the recurrence relation. We treat the preceding cases separately. 

A. u = n + ^. The required coefficients are given by 

where for j = 1, 2 

'' T{2k + 2 + 2l) 

1=0 



/32.+i.o(^) = L^2z(^) r(2fc + 2-20 

«i ( \ Y-^i ( , T{2k + 2l + 3) 
/32.+i.i(^) = LVi(") r(2fc + l_20 



1=0 

fc+i 



«2 , X V^ w , ,r(2A; + 3 + 20 



i=0 

fe+l 



r(2/fc + 2Z + 4) 
as well as that for 6ij the usual Kronecker symbol: 



p2 f \ — \^ JJ f \^ v'^"' -r ^t -I- 'i; 

/^2/c+2,ii^) -Z.^2mWr(2A; + 3-20 



r^j / ^ (-l)^r(fc + 2n + j + 2) ^ ,• ■ 

+ (1 - '^io)[7fc+,+i,o('^) - 7i+,-+i,i('^)] + Wk + j + 2n + 3) - log 2) + rfc+,-+i(i/)] 
where now for j = 0, 1 



fc+i 

72fc+2 



° oH = E^2/M^|§^fz^W2^ + 2fe + 3) + ^(2fe-2Z + 3)) 
72ViH = E4+iH^|§^|^W2^ + 2fe + 4) + ^(2fe + 3-2/)) 
72Vi,o(^) = E 4(^) ri2fc + 2-^o (^^^^ + 2fc + 2) + V^(2fe - 2^ + 2)) 
72Vi,i(^) = E 4+iH rJ2fct2-^0 ^^^^^ ^ ^^^ ^ ^^ ^ "^^^^^ ^ ^ " ^^^^ 

oo 

^m(i^) = E 4+m(^)r(^ + i)r(/ + 2m + 1) 

1=0 

We win use the above in order to get 52n+2fc+2(^)'-S2n+2fc+2(^)' C'2n+2fc+i(^)' C'2n+2fc+i(z')- 
II. 2z> ^ Z+. Let N = \2v] then we get that : 

Cl^k+ji^) = 2W+,+i \[-^^{-^f{lN+k+j,o(^)-^N+k+3,i^^))+PN+k+j+i{v)] 
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where 



rJV+fc+j, 



iN+k+jfii'') = Yl "^2 



1=0 



T{2iy + k + 2l + j + 2) 
21 y") Y{2v + k + j-2l) 






1=0 

oo 



P7V+fc+j+i(i^) = 5^ 4+N+k+j+2{^)^i} + N-2V + l)T{l + 2v + N + 2k + 2j + l) 



=0 



The preceding formulae wUl provide C0^+2fc+2('^)' C'2i.+2fc+i('^)' C'2r^+2fc+2('^)'C'2r^+2fc+2('^)• 
Certain elementary estimates. An apphcation of the steepest descent provides the following elemen- 
tary estimate for the F-function, for 5 € (0, 2)^^ — ^'• 



4S- 



csix + lfe~^^'^^^ < r(x) < Cs{x + i)^e-(i-^')^ 
while the Euler-MacLaurin formula gives for the '(/', -(/''-functions and a > 0, A; G Z^: 

7r2i(a) Ktpik + a + l)- log(/c + 1) + — — — r < 7rii(a), 

b[k + a) 

nuia) < V'(fc + a + 1) - (, ^\^^^^ ^) + ^^^2 < -22(a) 



as well as that 



,, , 1 4a + 3 , , ,, , 1 2a + 3 
7rii(a) = V(a) H h —, — — rr^, 7r2i(a) = V'(a) H h 



a 6(a + l)2' ''' ' ^^^' a 6(a + 1)2 ' 

, . ,/.. 1 3a + 5 a + 6 

7ri2(a) = V (o) n + TT, — — r\^ + 



a2 6(a + l)3 8(a + l)' 

,, .,, . 1 3a + 5 a + 6 

7r22(a) =ip [a) j + 



a2 6(a + l)3 8(a + l)' 
Actually we have the elementary inequality for i > 0, 



e"^ -1 
e^^- 1 



min{a,b) < —^ < max{a,b) 



n-^ = ^^ I fe-^'- 



which is used after the 

r(s) Jo 
These imply that for f3s < 1: 



- < dtii^) < 



\2l+u+i; - 'V '- p^lfj ^^^2l+u+^ 



that allow us to obtain that for Qim) = 2m + ^2 + 3^ r]i{5) = ^^^^zIt=^ , mi^) = ^'^V-^^^'^ *^^ estimates 
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K2(<5,l/)e('°g4-2-'?2('5))'» ^^(j^j,)g(2+log4-^i(5))m 
5 < rm{y),pm{v) < -—5 



{l + mY+2 {l+my 
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At this point we remark that the preceding identities are expressed in this notation in the form 



vr 



-{-lf{lN+k+jfi{^) - lN+k+j,l{^)) + PN+k+j+l{v) = 



sin 2vK 
The determinants of the recurrence relation are given respectively by 

Afc(z^) = [Bl^_^2k+2i^)Cln+2k+l{^) - Bln+2k+2{T^)C2n+2k+l{^)]i.n + k+ -)vl, 
^kil') = [Cly+2k+2i.^)Cly+2k+li.^) - Clu+2k+2{^)C2v+2k+li.^)\vl- 

The identities in conjuction with the elementary inequality 

e(l — x^) < logx < £{x'^ — 1) 
allows to exlcude the annihilation of the determinants. 

4.2 Case II. The coefRcients of the expansion 

In the preceding section we assumed that vq ^ Q when the 'initial power ' a = — 1. However we mentioned 
that when the initial power a > — 1 then we should use the Neumann series around the Bessel operator 
Hi^. The preceding paragraph suggests that certain terms are excluded from the first order perturbational 
terms and hence signify that the potential is decomposed as: 

v{x) = x^'^vq + x^'^u{x) + V 

where u G C^(R+) and as x — )■ 0, for U2j = V2u+2j then u{x) ~ ^1^=1 U2jx'^^ ■ Expanding the trace we 
concentrate on the contribution of the first three terms, abbreviating R\ = R\{H^) to obtain mod{X~°°): 

Tr{Rx ■ vRx ■ v) = 2uQTr{Rxx^'' Rx ■ v) + 2uoTr{Rxx'^''Rxx'^'' ■ u) 
+ ulTriRxx"^" Rx ■ x"^") + 2Tr{Rx ■ vRx ■ x^^u) 
+ Tr{Rx ■ vRx ■ v) + Tr{Rx ■ x^^uRx ■ x'^'^u). 

However, there we should assume that uo = V2u 7^ 0; otherwise we have through the commutator identity 

x^^Rx = Rxx^^ + 2jRxD2j-3x^''^-^^Rx 

which by iteration allows us to exhaust all the x— powers at the cost of i?A-powers. At this point the 
scaling argument restricts the contribution to a given order only to that from the first term. Therefore, 
the traces are essentially comprised in the general form 

I{a) := Tr{R^ ■ x^^ -Ra-V) 
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where V stands for x'^'^u,v. The identity 

d^x'+'l.+iix) = x''+H,{x)x 

suggests further that 

d^x^-^\ll+^{x) + Ilix)) = 2(2z. + l)x^'+^ll{x) 

and an integration of the latter reduce the study to that of integrals of the form for ^ = - : 

v2 roo 



where 



i{a) = ^ r ^x'-^^M-)HO 

^ Jo X 



Zu + L t, q, 



having set that 

The Singular asymptotics lemma suggests that if 2i/ + a ^ Z u q, 7^ is a nonvanishing asymptotic 
coefficient of the "potential " V then 

D2u+a+2,oI{0) = C2u+a+2{'^)uoVa + P2u+a+2{v) 

while if 2z/ + a G Z then 

D2u+a+2,lli0) = B2v+a+2{v)uQVa +Q2v+a+2{v). 

For these we need the Mellin tranfrom of J-{^) which are calculated in the next paragraph since 

C2v+a+2{'I^) = U2v+a+2fi{^) = T{-2v - a - 2), 
-B2y+a+2(i^) = -021^+0,4-2,0-^ = -ReS5=_(2i/+Q,+2) (-^)- 

4.3 Bessel function formulae. 

We give in some detail the calculation of the Mellin transform of the functions: Bo,Bi. The calculations 
of the Mellin transforms are executed using the integral representation for j = 0,1 

Jo X -\- t; X 

and the classical Weber-Schaftheitlin integrao. The convolution formula allows thus to obtain for: 

^^■^") = /r-wi^ ,2/0 -^ • / F,{w)F2{w)F3{w)F,{w)dw 

setting 

n^n-^ + u + j) ^^ r(-f)r(v + ^ + j) 



^iM := ' / ; , , ., , F2{w) : 



n-i + i^ + j) ' '' '■ r(-f + i. + i) 
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These are condensed in one formula: 



..- r , v.^ . ^^^dx _ r(f + |)r(^f + . + j) 



/ x^ ''J^+j{x)J^{x)— = — ^ 

Jo ^ J- (.9 



X r(is + i)r(f + i 
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ly 



Next we study the meromorphic properties of the function Bf the resultant formula implies in the domain 
D'i = {(u,v,vS) G (R+)^/t(; < 1, wi^ > u;} that: 



where 



(r(j + ^ - ^.))^ 






w , ,_ "w^ 



This is written further through a Taylor development as 

I.(s) = V ^(i + ^ + ^) J ( ^ 

if we introduce the integral 

/^,(S) :=/■/■ ^!^[^-.-2/-l(l +^2)-.+ f)][^-s+2.+i-2Z(;^ ^^2)-.-i+f][^2.^-(2.+l) ^ ;^]_ 



This represents a holomorphic function in the left half plane Ks < positive on the negative real axis. 
Furthermore the substitution u^ = ^37- — 1, f^ = ;^ — 1 leads to the familiar type of integral u 



-I rl 



his) = \ll^ri-^^W^^{i-o^^^-^{i-vO- 



|+^_;+i_l 



\1+ ^ ^ 



(l_^).(l_^^)-+i 



The following facts allow us to conclude that there exists a > — 2 such that one of the coefficients 
B2u+a+2{'^) or C2u+a+2{i^) is non zero. These are indeed checked easily: 

1. Notice that \Ii{s)\ < ^\Io{s)\. 

2. For s G C, Ks < 0, / G N we observe that \Iji{s)\ < Ijo{^s); precisely 

1 Li(s) 1 



(1 + M)2'+2^+l - Ijois) - (1 + M)2'+l+^ 

3. Actually as / —t- 00 then 

ra^S-^i) 

0{r 



r(f + i + _r^n-l-^^ 



l\T{v+\-l)T{v + 2 + ^ + 1) 
The series converges absolutely for all s £ C. 



'^lis = Au + joTS = 6u + j — 1 this integral is computed explicitly. We omit the result since this case falls out of our 
goals. 
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4- The function I{s) has poles at the negative of an odd integer. 

In the terminology introduced already we have the following lemma that has appeared in various forms 
in the literature, [4], [6] , [7] , [12] , [3] as as well as Melrose refers to it as the push forward lemma. 

Theorem. (The singular asymptotics lemma) Let f{y,x) G r°°(R^). Let {S,S') be the asymptotic 
character of /. Suppose that / has compact x— support and 

ms']f{y,x)\<{xyf^-^'hM 

for some S = 5z > and h^ satisfying f^ hzij)^ < oo. Let 

Jq X X 

Then F € r°°(R+) and Si + 5*2 is an asymptotic character of F. The asymptotic coefficients of F are 
given by 

5(fc)-l I 5'(fc)-l I 

r=j ■'' r=j 

- J2 ^'~"~^^'"' DlrDl,_i_J{y,x%=y=o 

where Ukj is a linear functional on r°°(R+) defined as follows. Let S be an asymptotic character of /. 
Let I > S{k) and Dz = xdx — z. Let 



nif){x) = fix) - Y. Y. Dzj f{0)x^ log- 

?tiz<?lik,Zf^kjeZ+ 



•■^ X. 



We define 



1 f°° dr 



It can be readily checked that Ukj{f) is independent of / and therefore of S, for S large enough. 
The computations are in fact facilitated through Mellin transforms because of the formula, [7]: 

Uk,j{f) = fj{-k), Dkj = /_j_i(-fc) 

where fi{zQ) is the coefficient of the {z — zq)^ in the Laurent expansion of the meromorphic extension of 

the Mellin transform: 

r°° dr 

f{z) = / f{x)x^- 

(which is defined for Rez >> if / is of bounded support) around z = zq. 
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